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We investigate the modulational instability and time-domain dynamics of nonlinear magnetic 
metamaterials composed of coupled split-ring resonators loaded by Kerr nonlinearity. Our results 
indicate that the recently proposed optical switching of local optical index based on uniform-response 
assumption seems fragile. We conceive two alternative schemes to utilize the valuable enhanced non- 
linearity, one is to focus on few-body systems and directly make use of the modulational instability 
(e.g., an optical comparator design), the other is to consider global switching arising from global 
feedbacks as in usual cases (e.g., Fabry-Perot cavities) rather than local-resonances-based switching 
of optical constants which may be destroyed by the discreteness and interactions as we show in this 
paper. We also try to provide comprehensive understanding of the relations between our results on 
discrete nonlinear metamterials and those on continuous metamaterials and natural materials in the 
literatures. 
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I. INTRODUCTION 

Linear metamaterials, ever since the theoretical 
works^^* and the experimental virorks^, have been high- 
lighted by a great number of achievements. Nonlinear 
metamaterials, hovifever, as Zheludev points out in 7], 
are sill "green apples" on the knowledge tree for meta- 
materials. 

In the discussions of nonlinear metamaterials, partic- 
ularly those composed of nonlinear split-ring resonators 
(SRR), it will be of great convenience to treat the meta- 
materials as continuous materials with controllable op- 
tical constants [e.g. ^{I) and n{I)], which allows us 
to employ existing concepts and methods in nonlinear 
optics. Several works since f§], including focus 
on this approach. They either directly assume intensity- 
dependent-index continuous media or assume uniform re- 
sponses of the discrete system before calculating optical 
constants. These approaches share a common implicit 
assumption, i.e., on a scale much smaller than the wave- 
length, a cluster of discrete inclusions in a metamaterial 
should uniformly (at least slowly-varyingly) respond to 
external radiation. This assumption justifies the con- 
tinuous descriptions which investigate the materials by 
partial differential equations. Particularly, some of them 
(0 and [iBl) propose applications for global-geometry- 
independent optical switching based on the bistable prop- 
erties of these local controllable optical constants. Note 
that these bistabilities^^i^ are obtained after assuming 
that all the meta-atoms respond uniformly in the mate- 
rials, at least in the equal phase planes, just as in the 
linear cases. 

However, we should note, that nonlinear metamateri- 
als (usually composed of nonlinear SRR) are intrinsically 
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discrete, and that they can even present strong local feed- 
back at the "single-atom" level, unlike the conventional 
nonlinear optical materials to which the feed-back is usu- 
ally introduced globally (e.g. by two mirrors in a Fabry- 
Perot cavity). Therefore, we expect to gain more knowl- 
edge if we treat the nonlinear metamaterials as discrete, 
interacting nonlinear lattices as they are. 

Some results have been obtained by treating the non- 
linear metamaterials as discrete nonlinear lattices, in- 
cluding What we should pay particular atten- 
tions is the series of work a^^d''''^^ !^^ , in which the authors 
show that the steady responses of dissipative SRR arrays 
may violate the translational symmetry and can be in the 
multibreather statesSS even when the arrays are driven by 
uniform radiation. The breather states in these works are 
theoretically constructed in the week-interaction regime 
and in the anti-continuous limilj^i in which the states are 
a priori nonuniform and may in practice be excited by 
introducing chirping radiatio n^^i^^ . 

In this paper, by investigating nonlinear effective cir- 
cuit models, we show that even when we start from 
the nearly-zero states under uniform external radiations, 
modulational instabilities will step in the bistable ranges 
proposed in [8|, i.e., the responses of the SRR array will 
spontaneously break the translational symmetry , and 
finally arrive at nonuniform states or even be chaotic, 
which will make the definition of optical constants in Q 
invalid at least in the bistable range of interest. These 
results reveal that the switching of field-dependent op- 
tical constants may not exist due to modulational in- 
stabilities. Despite the negative results, we propose the 
possibility of constructing all-optical comparators based 
on the spontaneous-symmetry-breaking phenomenon of 
two strongly coupled SRR. We also propose to utilize 
the valuable enhanced nonlinearity in SRR arrays, by in- 
troducing global feed-back mechanism as in usual cases. 
Finally, we comment on the relations of the results here 
based on our discrete models for metamaterials and those 
in the literatures. 
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II. THE THEORETICAL MODEL 

When a plane wave is normally incident into a three- 
dimensional nonlinear metamaterial composed of SRR, if 
we adopt the uniform assumption as in^, the SRR located 
on an equal-phase plane of the wave will experience uni- 
form external driving force (induced by other SRR and 
the external radiation) and respond with uniform ampli- 
tude and phase. We can thus simplify our considerations 
by focusing on a two-dimension (or even one-dimension) 
array of SRR radiated by uniform external fields. If we 
find that in some cases the uniform assumption is not 
valid even in the equal phase plane, we can undoubtedly 
conclude that such assumption is not valid in the three- 
dimension bulk materials. 

We model a rectangle lattice of SRR by LCR circuits 
coupled by mutual inductance (just as in|8|) with nearest- 
neighbor approximations and periodic boundary condi- 
tions. In terms of dimensionless quantities, the dynamic 
equation of the lattice under slowly- varying approxima- 
tion is 
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where u^j , g^j are the (dimensionless) local electro- 
motive force acting on the SRR at site (i, j) and the 
electric charge hold by its capacitor, respectively. The 
is the dimensionless frequency; t, the time; 7, the resis- 
tance; Kx.y, mutual inductances between nearest neigh- 
bors in the 'x, y' directions, respectively. Then the steady 
states should satisfy the equation 

Uij -I- K^n'^iqi^ij + qi+i,j) + Kyil'^{qij-i + qi.j+i) 
-{-n^~ + i-/n + l-\q\^)q,^, ^0 (2) 

More detailed definitions of the quantities are pre- 
sented in Appendix [Xj 



III. MODULATIONAL INSTABILITIES OF 
UNIFORM RESPONSES OF AN SRR ARRAY 

The modulational instabilities of continuous non- 
linear materials with negative refraction index have 
been studied in [TTj . in which temporal, spatial, and 
temporal-spatial modulational instabilities are all stud- 
ied. The spatial modulational instabilities of three- 
dimension discrete nonlinear SRR arrays induced by 
magneto-inductive waves have been studied in (21i] . in 
which it is shown that the unstable range may completely 
cover the decreasing range of the bistability, especially 
the two switching points. In this section, we in general 



employ the formalism in [21j but focus on the model ([T]) 
introduced in the previous section. 

Under uniform-response assumption as in Q, i.e., 
Wij = u and qij ~ q, the Eqs.© for steady states are 
simplified 



~i-in-l + \q\^)q = Q. (3) 



We can solve this equation and obtain the bistabilities 
of uniform responses as illustrated in Fig.(l). However, 
since mutual interactions in densely stacked metamate- 
rials are inevitably present, it is a nature question that 
whether possible excitations of magneto-inductive waves 
(spatial Fourier components of fluctuations) could induce 
modulational instabilities of these uniform states. By 
linear stability analysis (e.g. .[211 and [231] ). we can obtain 
the growth rate for each Fourier component of magneto- 
inductive waves Eq. (jB3p . If for all Fourier components, 
A(k) has a negative real part, the corresponding state is 
stable. Otherwise, unstable. 

After the linear stability analysis, the stable and insta- 
ble uniform responses are denoted by solid and dashed 
segments in Fig.lT]), respectively. Since we have chosen 
an equal phase plane to investigate, the mutual induc- 
tances between nearest neighbors should be positive in 
one direction and negative in the other direction. In all 
cases in Fig.([T|), we can clearly identify critical points as 
the boundaries between the stable and instable ranges. 
As the system is tuned from the stable ranges (across 
the critical points) to the instable ranges, we expect that 
it will spontaneously break the translational symmetry 
and actually be in nonuniform states. We can also ob- 
serve from Fig.([T]) and Eq. (jB3l) that in the low intensity 
limit, the modulational instability is absent for all cases 
implying that nonlinear metamaterials always allow uni- 
form description in perturbative regimes. 
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FIG. 1. The stability analysis of uniform response. For the 
curve with uniform bistability, = 0.835; for the other, Q. = 
0.850. Other parameters: = 0.3, Hiy = —0.1, 7 = 0.01. 
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IV. TIME-DOMAIN RESPONSE OF A 
TWO-DIMENSION SRR ARRAY 

A natural question arises that to what degree the 
actual response can deviate from the uniform assump- 
tion. In order to show an explicit picture of typical 
nonuniform responses when a system has evolved into 
the unstable range, we directly integrate Eqs.(Il]) using 
the Runge-Kutta method. In the recent papersi^iiiii^, 
single breather states (localized states which break the 
translational symmetry) are constructed from the anti- 
continuous limit in the weak-coupling regime. But these 
methods possibly can not be used in our problems since 
the coupling may be too strong for a single breather's 
existence. Here, we start from nearly-zero states with 
tiny random fluctuations and turn on the uniform radia- 
tion at the beginning. We then monitor the time-domain 
responses of two-dimension model (containing 10 x 10 
SRR) under different intensities of radiation. The results 
are shown in Fig. ([2]). The parameters are chosen as the 
case with uniform bistable response in Fig.(IT]). In that 
case, the critical radiation intensity u is about 0.0016 
which is the boundary between the stable range and un- 
stable range. In Fig.©, we test four radiation intensities 
|u| = 0.0014, 0.0016, 0.0018 and 0.020, of which the for- 
mer three are near the critical value while the fourth is 
deep in the unstable range. The initial response is nearly 
zero with each individual response randomly distribut- 
ing within [0,10~^]. For each test intensity, we evolve 
the system for 10^ periods and plot the time-domain re- 
sponses of the the SRR at two different sites, (1,1) and 
(5,5). The final states for each case are also plotted in 
3D figures. We observe that for the intensity below the 
critical value |m| = 0.0014, the response of the system 
quickly converge to a steady state which is uniform for 
all SRR, which means that the uniform assumption in Q 
is justified. However, when the intensity enters the unsta- 
ble range, e.g., |u| = 0.0018, the system spontaneously 
evolve into nonuniform states and cannot even reach a 
steady state (it seems like a limit circle state with peri- 
odic oscillation) . For the case very near to the boundary 
|u| — 0.0016, the situation is similar to \u\ = 0.0018 
with the translational symmetry also drastically broken. 
When we come to the case in the deep unstable range 
|u| = 0.020, it is even worse is that the time-domain 
responses quickly become chaotic. A common charac- 
teristic for the cases with nonuniform responses is that 
the actual states are highly disordered with individual 
responses being randomly positive and negative, which 
can hardly support regular wave propagation. 



V. OPTICAL COMPARATORS BASED ON THE 
SPONTANEOUS SYMMETRY BREAKING 

We can see from the above two sections that the re- 
sponses of nonlinear metamaterials in the unstable ranges 
can be rather complex. If we want to make use of these 
nonuniform states, apart from directly investigating the 
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FIG. 2. The time-domain dynamics for different radiation 
intensities. For each intensity the final states after r = 10* 
time of evolving are plotted in 3D figures. 

details of the complexity (e.g., [l3 and 111]) in these 
ranges, we may also focus on few-body systems to avoid 
these complexities. 

The spontaneous symmetry breaking (bifurcation) and 
its sensitivity on small deviations indicate that this phe- 
nomenon can be utilized to construct comparators. In 
FiglSl we illustrate this possible application by plotting 
the evolution of the steady states (obtained from Eq.([2])) 
of two coupled SRR under the continuous tuning of the 
overall intensities of the two electromotive forces which 
deviate from each other by a small amount (10"'^). Al- 
though in the initial state the responses of the two SRR 
only have a tiny difference, the resulting states in the 
symmetry-breaking range can have a large contrast ra- 
tio. This property can be used to compare which signal 
(in our case the electromotive force) has a larger initial 
amplitude. Due to the small size of SRR compared to the 
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working wavelength, the two SRR should be excited and 
detected via near-field ways (e.g. [13 )■ Such a device, 
if the two SRR can be fabricated as identical as possi- 
ble, can be used to realize ultra-sensitive and possibly 
utlra-compact all-optical comparators in applied areas. 
We should stress that such applications do not have to 
rely on SRR, and they may be even realized in any other 
devices if two identical nonlinear resonators are strongly 
coupled with each other, resulting possible spontaneous 
symmetry breaking (e.g., [28l - l30| ). These devices do not 
involve any assistance from electronic devices and thus 
may be of interest in the fields such as all-optical logic 
gates, all-optical signal processing and photonic comput- 
ers (e.g., |3ll - [33l |). One remarkable feature of our pro- 
posal is that we do not need high quality factor as other 
applications for optical signal processing2ii2iri21, since we 
only assume Q/ = 10 in our calculations. We also do not 
need uniform bistability. 
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FIG. 3. a): The principle of the proposed comparator made of 
two strongly coupled identical SRR. The solid (dashed) line is 
for ui — U2 — 0.001(— 0.001), respectively, b): The deviation 
from uniform states as functions of overall radiation intensity. 
The stability of the uniform states can be guaranteed because 
of the continuation procedure. Parameters: Q = 1.3, 7 — 0.1, 
K, = 0.6. 



VI. SETTING METAMATERIALS INTO 
FABRY-PEROT (FP) CAVITIES 

It is fascinating to see that due to field enhancement at 
the "single-atom" level, the nonlinearity of a single SRR 
can be strongly enhanced, and bistability can even occur 
for a single SRR at relatively low intensities, which cer- 
tainly is an advantage of nonlinear metamaterials com- 
pared to natural materials. If SRR are simply stacked 
together, however, problems will arise as we have stud- 
ied. One way out to utilize the enhanced nonlinearity 
is to focus on few-body devices as we studied in the 
previous section, i.e, directly making use of the multi- 
stabilities and modulational instabilities. Another is to 
operate the stacked devices under the low-intensity ra- 
diation and introduce global feed-back mechanism, e.g., 
setting them (e.g., the structure in 15]) in Fabry- Perot 
cavities with optical lengths much larger than the wave- 
length. This approach can not only effectively employ 
the enhanced nonlinearity originating from the sing-atom 



level, but also avoid the modulational instabilities suf- 
fered by the devices at higher intensities. That is be- 
cause the FP-bistability requires a jump of the refraction 
index only of order X/Lpp (compared to the negative- 
positive jump proposed in ^ and [isj ) , where A is the 
wavelength and Lpp is the length of the cavity. Detailed 
estimations reveal that due to the field enhancement in 
the narrow gap, the required radiation intensity inside 
the FP cavity would be smaller by a factor of about 
10^^ than that required if the FP is filled with the same 
kind of pure Kerr materials (see Appendix [C|). This SRR- 
enhanccd nonlinearity may be even more powerful and 
useful when combined with other nonlinearity-enhancing 
mechanisms, e.g., [ssl and Isoj . 



VII. DISCUSSIONS 

We want to comment on the relations between non- 
linear metamaterials composed of SRR and natural non- 
linear materials composed of natural atoms. The usu- 
ally mentioned bistabilities in natural materials are in- 
duced by a global feed-back mechanism and work at a 
perturbative regime, i.e., the nonlinearity is restricted 
to the third-order. In other words, usual optical bista- 
bilities are actually the nonlinear properties of certain 
global optical modes. No direct applications related to 
bistabilities in natural materials have been focusing on 
the non-perturbative regime in which the nonlinearity is 
very large at single-atom level or equivalently the radia- 
tion field intensity is comparable to the intrinsical field 
intensities in the atom level, since the responses may be 
very complex and very hard to control. Similarly, the 
proposals on controlling the bistabilities of the uniform 
optical constants^ii^ are just examples of working in the 
non-perturbative regime of the effective index of the bulk 
metamaterials, so it is not at all surprising to see the 
unpleasant instability problems. Actually, the uniform 
bistabilities in and fl5| are properties of local reso- 
nance responses rather global modes in usual case. Our 
proposal to introduce global feed-back by setting an SRR 
array in an FP cavity is just an attempt to avoid these 
non-perturbative regimes and try to seek a perturbative 
regime (using optical modes) to utilize the SRR-enhanced 
nonlinearity. We illustrate the difference between bista- 
bilities of local responses and global modes in FigH] 

We also want to stress on the relations between the 
spatial modulational instabilities discussed in this pa- 
per and those for continuous nonlinear media. Contin- 
uous descriptions assume homogenization of a material 
and employ partial differential equations. The modu- 
lational instabilities and related solitons in continuous 
descriptions (e.g., [HI), if indeed valid, are actually trig- 
gered by certain exploding long-wavelength (compared 
with the lattice constant) magneto-inductive waves^ 
(or equivalents in other discrete systems) and at the 
same time preserved by the decaying short-wavelength 
magneto-inductive waves in the discrete descriptions. 
However, our results in Fig.l and Fig. 2 show no indica- 
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tion of such long- wavelength-only modulational instabili- 
ties even in the low intensity range. All the modulational- 
instability-induced nonuniform states in Fig. 2 have short- 
wavelength components. This seems like a paradox. Re- 
member that the self-focusing-induced spatial modula- 
tional instabilities originate from the propagations of 
waves. As we state in the theoretical model section, we 
have only analyzed the modulational instabilities of an 
equiphase surface and have not involved the stabilities 
of propagations of waves, and therefore our results do 
not imply the paradox. These arguments implies that at 
comparatively low intensities, the spatial modulational 
instabilities can still happen due to propagation effects. 
As an added note, the modulational instability discussed 
in this paper should be termed as subwavelength modu- 
lational instability as in a recent publication [41]. 

We should also note that in the perturbative range 
(or low intensity range) in Fig[Tl the stable uniform re- 
sponses enable the materials to be perturbatively char- 
acterized as nonlinear media with Kerr nonlinearity (also 
indicated in the Appendix[C]). Therefore, the works (men- 
tioned in the Introduction section) which directly assume 
Kerr nonlinearities together with negative refraction in- 
dexes could be regarded as theories describing the pertur- 
bative properties of real discrete nonlinear metamaterials 
(e.g., as in^). 

One may question our negative results by referring to 
the experimental work S '-- "'^^ . Although frequency shift 
is demonstrated (certainly not sufficient for the index 
switching), no strong evidences (e.g., hysteresis curve) for 
bistability-related index switching are discovered, how- 
ever. Moreover, the nonlinear metamaterials in [43| and 
are excited by monopole antennas, a nonuniform ex- 
periment condition way too far from explicitly demon- 
strating the uniform bistability of optical constants. In 
fact, the two papers only claim that their metamateri- 
als are tunable and do not claim that the bistability is 
confirmed. Another question would arise from the LCR 
model we employ. Although this model is more or less 
simplified and is not general, it can capture the main 
features of many metamaterials, e.g., resonances, loss, 
mutual interactions, Kerr nonlinearity. Therefore, our 
results can give strong caveat for considering the appli- 
cations of the nonperturbative behaviors of effective op- 
tical constants (e.g., bistability of refraction index) when 
interaction-induced modulational instability in the lat- 
tice level is possibly present, since the effective optical 
constants in those cases may be ill defined. We also 
want to point out that the recent publication'*^ men- 
tioned above investigate another kind of nonlinear arrays. 
After making dipolar and slowly-varying approximation 
near a resonance, they derive almost the same formulism 
as in this paper and also obtain modulational instability 
in the uniform bistable range. 

We also feel it is necessary to clarify the differ- 
ence between the bistability of optical constants in the 
proposals^ii^ and the bistability induced by periodic 
structuresAS^'-^ in the literatures. The origin of the for- 
mer is the local feed-back at the single-atom level while 



that of the latter is just the periodicity, (termed as dis- 
tributed feed-back in the seminal paper—). The period- 
icity of the former has been homogenized, being replaced 
by a homogenized optical constant. What is more, the 
former is derived by local consistent equations of local 
optical constants and do not rely on the shape of the 
bulk material (explicitly shown in [8|; in [Tsj different 
shapes are used with the same interpretation), while the 
latter is derived from the global consistent equations ex- 
plicitly determined by the boundary or periodicity. One 
can also introduce such distributed or global feedbacks 
on a length scale larger than the wavelength to achieve 
bistability in metamaterials as we propose in the previ- 
ous section. These differences all arise from the working 
wavelengths of these phenomena, one is larger than the 
lattice constant while the other is comparable or smaller 
than that. Therefore, the historical success^i^ cannot 
be evidences for the theoretical proposals in^ii^. 

Finally, we should be reminded that if the whole array 
of SRR are rather sparsely arranged, i.e., the coupling be- 
tween neighbors can be neglected (gas phase), the local- 
bistability-based optical switching should also approxi- 
mately be valid although the two switching points are 
inevitably affected according to Ea. (jB3|) . But in those 
case the array is so sparse that it can hardly be treated 
as effective materials. The coupling parameters chosen 
in this paper is more typical in realistic cases (condensed 
phase) . 



VIII. SUMMARY 

By investigating the modulational instabilities and 
multistabilities, we find that some recently proposed op- 
tical switching based on the bistability arising from uni- 
form assumptions may not be valid. In order to get 
rid of this dilemma, we propose two approaches to uti- 
lize the valuable enhanced nonlinearity of metamaterials 
at single-atom level, one is to focus on the multistabili- 
ties of few-body problems (e.g., all-optical comparator) 
and the other is to operate them at low radiation in- 
tensities by introducing global feed-back mechanism to 
achieve global bistability effects (e.g., Fabry-Perot cav- 
ity) in perturbative regimes of optical constants. Finally, 
we briefly discuss the relations between our results and 
existing ones in the literatures. Particularly, we distin- 
guish the questioned bistability supported by local res- 
onance and the usual bistability of global optical modes 
(global resonances, or global feed-backs) and stress that 
the former may easily be destroyed by modulational in- 
stability induced by discreteness and mutual interaction 
in dense materials. 
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Appendix A: The Theoretical model 

The typical structure of SRR is shown in Pendry's 
seminal paper^. According to the effective LCR-circuit 
model of the SRR in the static limit, we can write down 
the dynamic equation of SRR lattice. 



-M. 



dt^ dt C 

<PQi+l,j d?Qi-l,j 



dt^ dt^ 

d^Qi.j+1 d?Qi,j-l 



dt^ 



dt^ 



(Al) 



where L, R and C is the effective self-inductance, re- 
sistance, and capacitance, respectively. The Qij is the 
charge hold by the capacitor at the site The M^^y 

are mutual inductances between nearest neighbors in the 
X (y) directions, respectively. Note that the nearest- 
neighbor approximation is employed. 

In our case where the split of the ring is filled by Kerr 
materials, the effective capacitance depends on the elec- 
tric field intensity. Eg, in the gap, i.e., C — Co(l -t- 
li^gP/liJcP). Here, the Cq is the zero-field capacitance 
and the \Ec\ is the characteristic field intensity. 

We also have the relation between Q and Eg in the 
quasi-steady states. 



Q 



CEgdg. 



(A2) 



where dg is the gap width of the capacitor. Then we 
arrive at the relation between Q and C, 



C = Co 1 



CUl\E,\^ 



Here since we usually work in the weak-nonlinearity 
regime of the Kerr material, the above relation can be 
approximated as. 



For the convenience of the discussion in Sec. (jVII) . we 
write down the nonlinear magnetic susceptibility of a sin- 
gle SRR. 



/ 1 
y oc — cx — . 

^ U Z 



(A6) 



where / is the current and Z is the complex impedance. 
With the linear part xi separated, it reads 



X = XI- 



-n^ + i-ffl + l - |e|- 



(A7) 



Appendix B: Linear stability analysis 



We study the linear stability of uniform states qij = 
q solved from the steady state equation ([3]). We first 
introduce tiny fiuctuations to the steady states, qij — 
q-\-Ai j, and derive the linear dynamic equations for j 
by use of the full dynamic equation. 



dAr, 



dr 



dA 



m-\-l.n > 



dr 



dAr, 



dA 



dr 



dr 



+ Arn,n+l) " {-^^ + ijQ + 1 - 2\q\^)Arn,n 
W^l^n (Bl) 

Since the lattice is uniformly periodic, the eigenmodes 
should be uniform plane waves. Thus let Am,n(k) = 

ae^k^m+tkyu ^ p^-ik^m-ikyu _ -y^g t\ien have 



d 
di 



a 



1 

irj 



A 

-B* 



B 

-A* 



a 

13* 



(B2) 



(A3) where A ~ 2kx^^ cos + 2Kyfl^ cos ky — (— fi^ -f- ijfl 
1 — 2|gp), B = q^ and irj ^ •y + 2iil{l + 2KxCOskx 



2Ky cos ky). The eigenvalues A of the matrix on the left are 
the growth rates of the magneto-inductive fluctuations 
with wave vector (kx, ky). 



and 



C « Co 1 + 



Cidl\E,\^ 



(A4) 



A(k) = - [~jn±{\qf - {2Kxn^coskx + 

2KyVl? COS ky + VL^ -1 + 2|gp)2)i/2] (B3) 



1 1 A \Q? 

C ~ Co ^ Cldl\E,\^ 



(A5) 



Combining the Eqs. (lAll IA5|) and make slowly- varying- 
amplitude approximation, we can obtain the dynamic 
equation ([T]) after defining the following dimensionless 
quantities, q = Cod^\E,\ ^ = I/VLCq, il = uj/ujq, 



7 = RCqujo, i^x,y = Mx^y/L, u = U/{\Ec\dg) and 



EJ\EA 



1- 



Appendix C: Estimations in Sec JVIl 

In the Sec lVn we propose bistability by setting an SRR 
array into an FP cavity. First of all, for the bistability to 
occur the required radiation intensity should reach a cer- 
tain scale to induce nonlinear deviation of the refraction 
index An ^ \/ Lpp, where A is the radiation wavelength 
in the metamaterial and L pp is the length of the FP cav- 
ity. This is usually a small factor provided that Lpp ^ A, 



7 



which will efTectively reduced the required radiation in- 
tensity and thus help to avoid the modulational instabil- 
ities at higher radiation intensities. In this regime, the 
nonlinearity of the metamaterial can be treated pertur- 
batively. 

Thus we need that the radiation propagate at least 
Lpp distance before it significantly attenuate. To this 
end, we should work at the frequencies where the meta- 
material is transparent, i.e., far from the absorption 
(1^1 = 10/Q/, where the detuning 5 — [uj — ujq)/ojq), 
and therefore the lossless approximation i? = can be 
justified. According to the results in the Theoretical 
Model section, the effective magnetization coefficient of 
the metamaterials can be approximately written in this 
case as x(|e|) = Xi[^ ~ |6p/(2(5)]. To express x in term 
of the radiation intensity we should relate the di- 

mensionless field intensity |e| in the gap to the radiation 
intensity \E\. 

The total electromotive force accumulated along the 
ring is 



y (ik X E) • dS| 



27r2r2 , 



A 



lEl 



(CI) 



with r being the radius of the SRR. Then in the lossless 
limit, the voltage V across the gap will be 



V 



U 



iujL + 1/ (iuiC) iuiC 
U 



U 



(C2) 



Now, considering the relations Eg ~ V/dg and e 
Eg/\Ec\, we have by relating the above results. 



\6\Xdg \E, 



Then x(l^^l) = Xi[^ ~ ^ib^lk;^]' which indicates that 
self-focusing and self-defocusing nonlinearity are both 
possible in the perturbative regime depending on the sign 
of the detuning S and xi ■ 

Now, we can estimate the radiation intensity required 
to induce the bistability of an FP cavity. 



4 4 



\E\' 



A 



XI 



(for the bistability without FP cavity, the right-hand side 
wiU be of 0(l).)then 



\E\V\Erf 



A m'^^di 

Lpp TT^r^x; 



Compared with an FP cavity filled with the same kind 
of pure Kerr material, the required intensity inside the 
cavity is reduced by a factor 



(C3) 



TT^r^Xl' 

For the structure proposed inii, if the Q / can reach 100 
(see Fig. 3a of^-), we can estimate S = 0.1. The expression 
for electromotive force (jCip and the effective-x descrip- 
tion should be still roughly valid for a wavelength-size ra- 
tio X/r as small as 5. The gap-size ratio dg/r is a tunable 
parameter depending on the fabrication technique which 
we assume can reach as small as 1/10. We also assume 
that the linear response (depending on the detuning and 
density of SRR) is not too small, i.e., xi ^ C'(l)- Then 
the reduction factor could be estimated to be about 10"*^. 

The role of the resonance can be assessed, from the 
above estimations, to be positive because it can ren- 
der large field enhancement [6 should better be small in 
Ea. (IC2p ]: negative because it inevitably bring in energy 
loss that may make the material opaque {6 should not be 
too small). 
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